Andreev reflection from non-centrosymmetric superconductors and Majorana bound 

state generation in half-metallic ferromagnets 
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We study Andreev reflection at an interface between a half metal and a superconductor with 
spin-orbit interaction. While the absence of minority carriers in the half metal makes singlet An- 
dreev reflection impossible, the spin-orbit interaction gives rise to triplet Andreev reflection, i.e., 
the reflection of a majority electron into a majority hole or vice versa. As an application of our cal- 
culation, we consider a thin half metal film or wire laterally attached to a superconducting contact. 
If the half metal is disorder free, an excitation gap is opened that is proportional to the spin-orbit 
interaction strength in the superconductor. For electrons with energy below this gap a lateral half- 
metal-superconductor contact becomes a perfect triplet Andreev reflector. We show that the system 
supports localized Majorana end states in this limit. 

PACS numbers: 74.45. +c,74.78.Na, 75.70.Cn, 75.70.Tj 
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I. INTRODUCTION 

Heterosystems with adjacent superconducting and fer- 
romagnetic phases may show unconventional spin-triplet 
superconducting proximity effects even if the supercon- 
ductor is of the conventional s-wave spin-singlet type.— 
Triplet correlations, even if they are weak, are important 
in ferromagnets, where the standard spin-singlet prox- 
imity effect is short-ranged as a result of the exchange 
splitting d In half- metals singlet pairings are ruled out, 
since in a half metal one spin species has zero density 
of states at the Fermi level, so that the triplet version 
is the only possible form of the superconducting proxim- 
ity effect. Microscopically, the superconducting proxim- 
ity effect is mediated by Andreev reflection, the phase- 
coherent reflection of an electron into a hole or vice versa 
at the superconductor interfaced Triplet superconduct- 
ing correlations then require a form of Andreev reflection 
that includes a spin flip.— £ 

The triplet proximity effect has been considered first 
in ferromagnets (with a finite minority spin popula- 
tion), where the observation of long-range superconduct- 
ing correlation effects^— has been shown to be consistent 
with the existence of induced triplet correlations in the 
ferromagnetic A number of mechanisms that give rise 
to the spin-flip Andreev reflection required for the triplet 
correlations, such as magnetic domain walls^ spin-orbit 
interaction^ and unconventional pairing correlations, 1- 
have been studied in hybrid ferromagnet-superconductor 
systems. 

Long-range proximity effects have also been observed 
in half-metallic Cr0 2 -NbTiNH and Cr0 2 -MoGei* het- 
erostructurcs. Because of the absence of minority car- 
riers in the half metal it is concluded that the observed 
superconducting correlations must be of triplet typej 13 ' 15 
However, in half metals the conditions for spin-flip An- 
dreev reflection are more restrictive than in a ferromag- 
net. In particular, electron-hole symmetry and cur- 
rent conservation pose stronger restrictions on candi- 
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FIG. 1: Top panel (a): Interface between a half metal (H) 
and a superconductor (S). For a superconductor with spin- 
orbit scattering, normal reflection as well as triplet Andreev 
reflection — the reflection of a majority electron (e) into a ma- 
jority hole (h) — take place at the HS interface. Bottom panel 
(b): A lateral contact between a half metallic film or wire of 
thickness d and a superconductor. For a lateral contact, an ef- 
fective Andreev reflection amplitude rjf for electrons moving 
towards the contact region can be defined, which describes the 
combined effect of multiple reflections at the H-S interface. 



date mechanisms for spin-flip Andreev reflections than 
in ferromagnetsi 16 ' 17 If the interface is symmetric with 
respect to reflection in the surface normal, electron-hole 
symmetry and unitarity require the Andreev reflection 
amplitude r\ le (e) to vanish at the Fermi energy e — 
Thus, mechanisms that give rise to triplet Andreev reflec- 
tion must either break electron-hole symmetry, i.e. take 
place away from the Fermi energy, or orbital symmetries. 
Mechanisms of the latter type are magnetization gradi- 
ents in the half metal ) 17 i 18 or impurity scattering^ 

In this article, we study spin-orbit interaction (SOI) 
in S as a possible mechanism giving rise to spin-flip An- 
dreev reflection in a half metal (H ). The presence of SOI 
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is contingent on the breaking of inversion symmetry. Ex- 
amples for systems with SOI in the normal state are sur- 
face states in Au,— semiconductor heterostructures and 
2D electron gases in quantum wells with partially tun- 
able SOIsi2i~— Crystalline superconductors such as the 
non-centrosymmetric cuprates and heavy-fermion com- 
pounds such as CePt3S*2i and others,— show SOI due 
to the absence of an inversion center in their crystal 
structure. These materials have received intense inter- 
est because they display unconventional superconduct- 
ing (helical^r— ) phases with mixed singlet and triplet 
pairing correlations ) 25 ' 29 " — magnetoelectric effects ) 32 ' 33 
and an anisotropic spin susceptibility 34 ' 35 But also cen- 
trosymmetric superconducting materials show SOI in 
surface or interfacial layers with inversion asymmetry. 
For instance, the breaking of inversion symmetry at a 
plane interface due to a change in the chemical potential 
gives rise to the Rashba SOI . 21 ' 34 

In our calculation we assume a general model for the 
SOI that is linear in momentum and includes the Rashba 
SOI as a special case. In the first part of this article, we 
derive expressions for the electron-to-hole and hole-to- 
electron Andrcev reflection amplitudes and r c h in a 
model Hamiltonian with parabolic dispersion in the half 
metal and in (the normal state of) the superconductor 
and to first order in the SOI. For our model Hamiltonian, 
which includes effects of a Fermi-velocity mismatch and 
a tunnel barrier at the half-metal-superconductor inter- 
face, we find that the Andrcev reflection is such that the 
induced superconducting correlations in the vicinity of 
the superconducting interface arc of even-frequency and 
complex p-wavc type. This is a significant difference with 
other possible sources of triplet superconducting corre- 
lations in half metals, such nonuniform magneti- 
zation direction in the half metal, which also allow for 
odd-frequency s-wave proximity effects. 

Because Andrcev reflection relies on the presence of 
SOI in the superconductor, the Andreev reflection prob- 
ability | rho | 2 may be small depending on the strength of 
the SOI such that the induced superconducting correla- 
tions become weak. However, a fully developed proximity 
effect 3 ^ can be achieved in a geometry in which multiple 
Andrcev reflections occur. Examples of such geometries 
are a half-metallic film or wire in contact to a super- 
conductor (shown schematically in Fig. [1] b) . The latter 
example is closely related to recent proposals for the real- 
ization of Majorana fermions in a solid-state settingj^Zr— 
which in turn play an important role in proposals for 
topological quantum computation . 37 ' 44 The second part 
of this article contains an investigation of multiple An- 
dreev reflections in the film- or wire geometry. We will 
show that the effect of multiple Andreev reflections can 
be combined into effective Andrcev reflection amplitudes 
r^f and r^f . These effective amplitudes may have unit 
magnitude if sufficient reflection events contribute coher- 
ently them. We show that in this case localized Majo- 
rana states can be formed at the ends of a half metal wire 
in contact to a superconductor with spin-orbit coupling. 



The investigation of the thin-film geom etry is motivated 
by the recent experiments of Rcfs. Il3lfl3 . 

In the present calculation we do not consider scattering 
from impurities in Hor S. This does not seriously affect 
the first part of our calculation, which addresses the An- 
dreev reflection amplitude for a single reflection off the 
HS interface, because the Andreev reflection amplitude 
is a local property of the interface. It is, however, a lim- 
itation for the second part of our calculation, since the 
proximity effect induced by the multiple Andreev reflec- 
tions is of p-wave type, which is suppressed by impurity 
scattering. Thus, the results derived in the second part 
of this article are valid only if the elastic mean free path 
in the half metallic film or wire is sufficiently large. 

The article is structured as follows. In Sec. Ill Al we in- 
troduce the model Hamiltonian of the HS heterostructure 
and describe the calculation of the Andreev reflection 
amplitudes for a single reflection off a HS interface. As 
an application of our calculation, we calculate the sub- 
gap conductance of a half-metal-superconductor inter- 
face and the Josephson current in a superconductor-half- 
metal-superconductor junction in Sec. IIII1 In Sec. II VI we 
consider the geometry in which a thin half-metallic film is 
brought into electrical contact with a superconductor and 
derive the effective Andreev reflection amplitude for 
this situation. Applications to the subgap conductance 
and Josephson effect in lateral HS and SHS junctions are 
then given in Sec. |Vj Finally, in Sec. IVII we consider a 
half-metallic wire placed in contact to a superconductor 
and relate our findings to predictions of the occurrence 
of Majorana fermions in such a system. 



II. INTERFACE BETWEEN HALF METAL AND 
SUPERCONDUCTOR 

A. Hamiltonian 

We consider the interface between a half metal (H ) and 
a superconductor (S) as shown in Fig. Q] a. Coordinates 
are chosen, such that the superconductor and the half 
metal occupy the half spaces z > and z < 0, respec- 
tively. Electron and hole excitations at excitation energy 
e near the interface are described by the Bogoliubov-de 
Gennes equation 



(1) 



where ^ = (u-f, u±, v-f, Uj,) is a 4-component wave func- 
tion with separate amplitudes for the particle/hole ex- 
citations («, v) in the spin up/down bands (f, i)- The 
Bogoliubov-de Gennes Hamiltonian H. has the general 
form 



n = 



Ho ia 2 A(r) 
- i( r 2 A(r)* -HZ 



(2) 



Here, the superconducting order parameter A(r) = 
Aoe^O(z), where Q(z) = 1 if z > and otherwise, 
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and the at are the Pauli matrices, i = 1,2,3. We take 
the normal-state Hamiltonian H to be of the form 

H = ^-- M*)£r + hw5(z) + H so , (3) 



<r=t4 



where m is the electron mass (taken to be the same in H 
and S), 



jU H (7 if z < 0, 
/is if z > 0, 



(4) 



with the potentials /iHf) Mh^, and /ig representing the 
combined effect of the chemical potential and band offsets 
for the majority and minority electrons in the half metal 
and for the superconductor, respectively, and where w 
sets the strength of a delta-function potential barrier at 
the interface. The operators 



A =5(1 + ^3), P i =\{l-az) 



(5) 



project onto the majority and minority components, re- 
spectively. (The magnetization direction in H is taken as 
the spin quantization axis, which need not coincide with 
the z-axis.) We will take the limit /xhj. —00, such that 
only the majority spin band is present in H. We further 
write 



2m 



2 1,2 



Ms 



P.S 



2m 



(6) 



where £;f,h and fcp.s are the Fermi wavenumbers in H 
and S, respectively, and 



mi 



2m 



(7) 



where k is the minority wavefunction decay rate. The 
Fermi velocities arc defined as 



i> FjH = ^fcp iH /m, v s ,n = hks,n/m. 



(8) 



The step function model for the superconducting order 
parameter A(r) is justified for s-wave superconductors if 
the coupling to the half metal takes place via a tunnel 
barrier with transparency r <C lr^ which corresponds to 
the requirement that \w\ ^> fp,H, wf,s- 

The operator i?so represents the effect of spin-orbit 
coupling. We consider the case that Hso is linear in the 
momentum p and that Hso is nonzero in S only— 



H so = ~[pe(z) + o(z)p]-J2nm, (9) 

where we denote = (Qi )X , &i, v , Qi yZ ) T ■ Such a SOI 
may originate from the breaking of inversion symmetry 
by the crystal structure of S or due to an inversion asym- 
metry of the HS heterosystem. We assume that the spin- 
orbit interaction is weak, ?L|0(1cf,s) I ^ v f,h, vs,n, so 



that it can be captured by treating -ffso to first order in 
perturbation theory. 

In addition to the spin-singlet order parameter con- 
tained in Eq. ©, the presence of SOI generally allows for 
a triplet contribution to the order parameter, which is of 
the form A(p) = ^i(p) cr «- Because of the Pauli 

principle, these triplet components are odd in momen- 
tum, Aj(— p) = — Aj(p). Such a triplet contribution is 
absent if the pairing interaction is isotropic^ but it may 
be present if the pairing interaction is anisotropic^ In 
appendix [B] we include triplet pairings in the model and 
give the results for the Andreev reflection amplitudes. 
The spin-orbit interaction does not lead to a modification 
of the magnitude of the spin-singlet order parameter to 
first order in Hso- 



B. Andreev reflection amplitudes 

We now calculate the Andreev reflection amplitudes 
for of the interface between H and S using the Blonder- 
Tinkham-Klapwijk formalism^ At the HS interface 
triplet Andreev reflection occurs because quasiparticlcs 
incident on the interface from H penetrate the supercon- 
ductor over a finite length before being reflected. Due to 
the SOI, spin is not a good quantum number in S, which 
makes spin-flip reflection possible. The Andreev reflec- 
tion amplitudes are found by matching eigenfunctions of 
H in H and S to linear order in the SOI. (An alternative 
method, using perturbation theory in the SOI Hamilto- 
nian, will be described at the end of this subsection.) 
The matching conditions, continuity and conservation of 
particle flux, hold for plane-wave eigenstates in the im- 
mediate proximity of the interface on length scales of the 
Fermi wavelength. Thus, the S-matrix of the interface is 
a local property and will not be changed by weak disor- 
der. 

Starting point of the matching procedure are expres- 
sions for the general solutions of the Bogoliubov-dc 
Gennes equation in Hand S, near the HS interface. Be- 
cause of translation symmetry along the interface, we can 
consider plane- wave solutions with wavenumbers k x and 
k y in the x and y directions parallel to the interface. In H, 
one then finds six linearly independent solutions, which 
we label * e ,t,±, *h,t,±, *c4> an d 



p ±ik z (-\-s)z-\-ikr E x-\-ik y y 

*c, t ,±(r) = == (1, 0, 0, 0) T , 



*h,t,±(r) 



e ^fik z (-£)z+ik x x+ik y y 



(0,0, l.or 



*h,4.(r) = e K * { - e)z+ik * x+ik i> y (0, 0, 0, 1) T , 



(10) 

(11) 

(12) 
(13) 



where k z {e) and n z (e) are the positive solutions of 
k z {e) 2 = fc|.H- fc(|+2me/a 2 , 
n z (e) 2 = k 2 + fen - 2m,e/h 2 , (14) 
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and 

k|| =(k X ,ky,Of, (15) 

is the momentum parallel to the interface and 



v z (e) = Kk z (e)/m. 



(16) 



The states labeled with + and — are majority states 
moving towards or away from the interface, respectively. 
They are normalized to unit flux. The states labeled with 
J, are minority states that decay into the half metal. They 
appear in intermediate stages of the calculation only and 
their normalization is not important. 

Only the spin-orbit interaction terms proportional to 
Oi and Q2 give rise to spin flips in the superconduc- 
tor. For a calculation of Andreev reflection amplitudes 
linear in the SOI, it is then sufficient to set O3 = 0, 
which significantly simplifies the form of the solutions 
of the Bogoliubov-de Gennes equation in S. In S, one 
then finds four linearly independent solutions ^f a> t of the 
Bogoliubov-de Gennes equation, t, s = ±1, which read 



/ 



*.,*(*) 



\ 



-i(j>— isrj-\-ij Stt 



V e 



-i(p—%srj 



Jqs^z+ik^x+ikyy 



where 



i] = arccos(e/A ), 



. (17) 



(18) 



Qs,t = t\j s — fen + 2itmJ Aq — e 2 — msif2^ t , (19) 



and 7 Sj t and f2 Sj t are defined such that 
(n x + Oh) ■ q = ^s,tqe t% 



(20) 



with q = (k x ,k y ,q St t) T - 

A complete solution \&(r) of the Bogoliubov-de Gennes 
equation consists of a linear combination of the six spe- 
cial solutions (fTU |) - (fT5|) in H for z < and a linear com- 
bination of the four special solutions (jTTJ) in S, with the 
boundary conditions 



*(r) 
d#(r) 



z=0~ 



z=0+ 



(21) 



dz 



<9#(r) 

dz 



=0+ 
3 



I /^<> ;: o\ + ^ ] *(i 



z=0+ 



(22) 



at the interface z = 0. Since ^(r) is a four-component 
spinor, the boundary conditions provide eight linear rela- 
tions between the coefficients of the ten basis functions. 
The Andreev amplitude rj lc is then defined as the coef- 
ficient of ^h,t,- if t ne coefficients of the two incoming 



wave solutions ^ e ,t,+ an d ^h,f,+ are chosen to be 1 and 
0, respectively. Analogously, the Andreev amplitude r e h 
is defined as the coefficient of ^e.t,- if the coefficients of 
^e.t.-f- an( f ^h,t,+ are chosen to be and 1, respectively. 
To lowest order in the SOI and to lowest order in the 
normal-state transmission t(9) of the HS interface, we 
find that 

^mT(g) e -^(ft 1+ z0 2 ).k||Ao 
rhe(k||,e)=-^ . nrn — — ( 23 ) 



F,S 



fc F>H sin^ 



and 



eh(k||,e) = r he (-k||,-e)* 

-miT(fl)e'* (n t - m 2 ) ■ k|| A 



2(fc 



F , s -fc2 S in 2 0) yA2— 72 



(24) 



Here 8 = arcsin(|k|| |/Af,h) is the angle between the in- 
cident momentum and the interface normal, see Fig. [TJ 
and 



, m2 W F,H C0S 0(«F,S- 
r W = 3 



vl „ sin 2 1 



1 



w 



(25) 



Equation (j2"5)l has been simplified using the "Andreev ap- 
proximation" , which amounts to neglecting corrections of 
order C(e/£ , F l s, Ao/-Ef,s)- (This approximation is uni- 
formly valid for all angles if &f,s > &f,h- If ^f.s < &f,h 
there is a small range of angles for which the approxima- 
tion fails.) 

The divergence for e -> ±A in Eqs. ([221) a nd (JMJ) is a 
consequence of the expansion in the normal-state trans- 
mission coefficient r of the HS interface and has to be 
cut off for 1 — (e/A) 2 < r 2 . This means that the imme- 
diate vicinity of ±A has to be excluded from the region 
of validity of Eqs. (|23| and (|24| . so that these equations 
are valid for 1 — (e/A) 2 ^> r 2 only. The same condi- 
tion will be required for the validity of Eq. (|26p below 
and for expressions that are derived from these equa- 
tions. (We note that similar restrictions also apply to an 
expansion in the transmission coefficient for a normal- 
metal-superconductor interface, see, e.g., Ref. l48l.) 

For completeness, we also give the results for the nor- 
mal reflection amplitudes r ee and r^h of the HS interface 
consistent with the assumptions of our calculation 



:(k|i,e) = -l 




(26) 



r h h(k||,e) = r* c (-k| h -e 



(27) 



Alternatively, Eqs. (|2"5)l and (|24p can be obtained from 
a calculation of the first-order perturbation theory cor- 
rection to the scattering matrix of the HS interface with- 
out spin-orbit interaction. This calculation is outlined in 
the appendix. (See Ref. [Tj] for more details.) 
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Equations (|23|) and ([24)1 are the two central results of 
the first part of this article. Although the Andreev re- 
flection amplitudes have been derived for a specific model 
Hamiltonian and in the limit of a tunneling interface, we 
believe that the symmetry properties of The and r e h — Hie 
and r e h are odd in k|| and even me — persist in a more 
general calculation, as long as the SOI is linear in mo- 
mentum. We have verified this statement for the cases 
that a finite minority-wavefunction decay rate k is in- 
cluded in the calculation, that the spin-orbit interaction 
extends only a finite distance into the superconductor, 
and that higher-order terms in the interface transmission 
t(9) are included. (See App. [2] for details.) 

The antisymmetry of rh c and function of ky 

implies that the Andreev reflection amplitudes rh and 
r c h contain only four elements fii^, fl2, x , and O2, y 

of the spin-orbit coupling matrix. We had already dis- 
cussed, that the three elements ^3^, &3,y, and fls :Z that 
describe the coupling between the spin component par- 
allel to the magnetization direction and the orbital mo- 
tion of the electrons do not give rise to spin flips and, 
hence, do not contribute to the Andreev reflection am- 
plitude. Equations (j2"3")l and (f!H|) show that the same is 
true for the elements £li tZ i ^2,z, and ^3 ^ of the spin- 
orbit coupling matrix that couple the electron spin to 
the orbital motion perpendicular to the interface and, 
thus, provide a spin- flip mechanism that is symmetric in 
k|| . For zero excitation energy e, this observation can be 
understood from the general symmetry considerations of 
Ref. [I?], which state that rh e (k||,0) = if ri lc is a sym- 
metric function of k|| . That this remains true for nonzero 
e is special to the case of spin-orbit coupling as a source 
of spin- flip scattering and requires the explicit calculation 
of this section. 

There is a direct relation between the Andreev reflec- 
tion amplitudes and r h and the anomalous Green 
function /(k, iw),— 



A. Subgap conductance 

We assume the interface to have lateral dimensions 
W x x W y and impose periodic boundary conditions in 
these directions. This leads to a quantization of the 
transverse modes with wave numbers fc„ x = 2im x /W x 
and k ny — 2Tm y /W y , n x and n y integer. At zero tem- 
perature, the differential conductance G = dl/dV can be 
calculated in terms of the Andreev reflection amplitudes 
The- Replacing the summation over modes by an integral, 
we findifi^a 



G(V) 



2e< 



Tr|r he (k||,e^)| 



(29) 



where Tr{...} = W x W y /(4ir 2 ) <2k„ {. . .} is the 

trace over transverse modes k|| . The factor two is due to 
the doubling of the transferred charge by conversion of 
an electron into a hole upon Andreev reflection. 

Substituting Eq. (|23l) for the Andreev reflection ampli- 
tude and performing the integrations over k x and k yi 
we then find 



G(V) 



3N (r 2 ) £ 



A 2 
^0 



h 8 



A 2 



(eV) 2 



(30) 



F,H 



where N = fc§ H W x W y /47T is the number of propagating 
channels in H and 

tt/2 

d6 . . . sin 9 cos 9 (31) 



denotes an average over the angle of incidence 9. 



f(k,iu) oc 



Q(-k z )reh(k\\,ioj) if to > 0, , * 
-e(A,)rhe(k||, -««;)• ifw<0, [M> 



up to a prcfactor that is not important for the identi- 
fication of the symmetries of /. Since fhe(k|i , — ilo)* = 
— r h(k|| , iuj) in the present case, see Eq. (|2~i)) , one con- 
cludes that the induced superconducting correlations in 
H are odd in momentum (i.e., predominantly of (com- 
plex) p-wave type) and even in frequency. 



III. APPLICATIONS: SUBGAP 
CONDUCTANCE AND JOSEPHSON CURRENT 



As an application, we now calculate the subgap con- 
ductance of an HS junction and the Josephson current 
of a supcrconductor-half-mctal-superconductor (SHS) 
junction. 



B. Josephson current 

For the Josephson effect, we consider a half- metallic 
junction of length Lj separating two superconducting 
contacts. Again, we take the junction to have lat- 
eral dimensions W x x W y and impose periodic bound- 
ary conditions in the x and y directions. Taking periodic 
boundary conditions is justified if the lateral dimensions 
W XiV Lj, see Fig. [3] a. We further take both HS inter- 
faces to have the same normal-state transmission t(9), 
take the same spin-orbit interaction Hamiltonian Hso in 
both superconductors, and neglect impurity scattering in 
the half-metallic junction. 

The Josephson current can be found from the den- 
sity of states in the junction which, in turn, may be 
expressed in terms of the scattering matrices of the HS 
interfaces (see Ref. [50II5H for details) . In the absence of 
impurity scattering, the contributions to the Josephson 
current from different transverse wavevectors k|| add up, 
and one finds that the Josephson current at temperature 
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FIG. 2: Serial (top, a) and a lateral (bottom, b) 
superconductor-half-metal-superconductor (SHS) junction. 



T is given by the expression 
2ek B T d 



I 



E 

"i~ n=0 

(Trjlndet [l - K(k h iuj n )K' (k h iuj n )] |,(32) 



where w„ = (2n + l)nkBT is the Matsubara frequency, 
1Z is a 2 x 2 reflection matrix containing reflection and 
transmission amplitudes for the first HS interface, 



e ifc i (icj„)Lj 

r ee (k||,iw„) r e h(k||,iw n ) 



' cc v (I J " w n / ' en v |p t "- w ri / \ /qq\ 

r hc (k||,iw„) r hh (k||,iw„) y ' 



with fe z (e) is given in Eq. (| 14[) and the reflection ampli- 
tudes given by Eqs. (|23j) - (j27)) above, and TZ'(k\\ , iuj n ) is 
a similar matrix describing Andrccv reflection at the sec- 
ond HS interface. Specifically, in the serial geometry in 
Fig. [2] a, 7£'(k||, iw n ) = 7£(k|| , io; n )|^ = o with the phase of 
the order parameter set to zero. 

Closed-form expressions for the Josephson current I 
can be obtained in limiting cases only. In the limit of 
"long" junction, Lj ^> £, where £ = Hvf^/Aq is the 
superconducting coherence length, and for high temper- 
atures UbTL-JTivy,^ 3> I one finds 



2e 2 

e7 =—j— 87rfcBTsin(0) 



(34) 



x Tr 



{|r h e(k |h ia;o)|| =0 e^^(o)} 



Performing the integration over the transverse momen- 
tum k|| for the parabolic dispersion of our model Hamil- 
tonian then gives the result 

el = 2wkBT 3 ahl4 ' G(0)f(2nk B TL i /hv F , R ), (35) 

where G(0) is the zero-bias conductance of a single HS 
interface, see Eq. (|30j) . and 

f(x) = e" x (6 - Wx - x 2 + x 3 ) + x 2 (x 2 - 12)Ei(-a?) 
KA8e~ x /x 2 forx>l, (36) 



with the exponential integral Ei(a;) = — dte~ t /t. 

In the opposite limit of zero temperature, the expres- 
sion for the Josephson current in a long junction (Lj ^> £) 
becomes 



2e sin t 



duo 



(37) 



x Tr 



f \n^l 

\ cosh - Re [r* e e? ik ,(p)L] + | rhe |2 C0S( 



Normal reflection with amplitude r co at the two HS in- 
terfaces gives rise to terms in Eq. (|37|) that oscillate 
with the junction length L y These oscillations disap- 
pear once the trace over transverse modes is taken, since 
k z (Q)Lj ^> 1 for long junctions. The remaining non- 
oscillatory contribution to the supercurrent / can then be 
calculated by taking the average I = (2t:)~ 1 dxl(x), 
where I(x) is obtained from Eq. (|3"T)) by the replacement 
2ik z (Q)Lj — > x- One thus obtains 



e sm m „ 
I = -^1^1^ log 



16 sin 2 ((j>/2) 



|r hc | 2 sin 2 i 



(38) 



The remaining trace over modes can be performed to log- 
arithmic accuracy by neglecting the dependence of the 
argument of the logarithm on k|| . This amounts to the 
replacement |r hc | 2 — > (|r hc | 2 )e = hG(0)/2e 2 N in the ar- 
gument of the logarithm. One then obtains 



el 



15L, 



G(0) sinolog 



32e 2 iVsin 2 (0/2) 



hG[0) sin 2 (j) 



(39) 



up to corrections of order G(0)Hvf / Lj, but without the 
large logarithm log(e 2 N/G(0)h)). The small corrections 
to the approximately sinusoidal phase dependence of the 
supercurrent in Eq. (|39[) originate from scattering pro- 
cesses with multiple normal reflections at the HS inter- 
faces. 



IV. LATERAL GEOMETRY 

An experimentally relevant situatio n 13 ' 14 is the lateral 
geometry where the superconducting contact is attached 
laterally to a thin H film. This situation is shown in 
Fig- lb- In comparison to the serial contact considered in 
the previous section, a lateral contact has a significantly 
larger contact area per unit cross section of H. Multiple 
reflections occur at the H-S interface, because quasipar- 
ticles are repeatedly reflected backwards from the lower 
film boundary towards the interface. In the absence of 
impurity scattering in the half-metallic film, the coherent 
addition of these multiple Andreev reflections leads to a 
significant enhancement of the Andreev reflection prob- 
ability for a quasiparticle incident on the lateral contact 
from the left (in Fig. QJ>), as we now show. 

We choose coordinates, such that the half metal oc- 
cupies the region between z = and z — —d and the 
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superconductor occupies the region x > 0, z > 0, see 
Fig. [TJd. We take periodic boundary conditions in the y 
direction, with system size W y . For the half metal, we 
take hard-wall (Dirichlet) boundary conditions at z = —d 
for all x, and at z = for x < 0. The thickness d of the 
half-metallic film is taken to be much smaller than the 
superconducting coherence length £ = hvp.n/ Aq. As be- 
fore, we take the HS interface to be a tunneling interface 
with a transmission probability t(9) <C 1. 

The goal of our calculation is to find the amplitude 
r^f that a right-moving electron-like quasiparticle ap- 
proaching the contact from the left is Andreev reflected 
into a left-moving hole-like quasiparticle, as well as the 
amplitude for the process that a hole-like quasipar- 
ticle is Andreev reflected as an electron-like quasiparti- 
cle. The calculation proceeds in three steps: First, we 
construct scattering states in the absence of spin-orbit 
interaction; Second, we account for the effect of spin- 
orbit interaction in a superconducting region of length 
d <C 8L <C hvp t ji/Ao using perturbation theory; Finally, 
we combine Andreev reflections from different segments 
and compute the Andreev reflection amplitudes r£f and 

r cS 
'eh ■ 



A. Scattering states in the absence of SOI 

Because of translation symmetry in the y direction, 
the scattering states can be chosen as plane waves in the 
y direction with wavenumber k y , which takes discrete 
values only because of the periodic boundary conditions 
in the y direction. We first construct scattering states 
for x < 0. There, because of the hard-wall boundary 
conditions at z = and z = —d, the z-dependence can 
be chosen proportional to sin(fc z z), where k z = nir/d, 
n = 1,2,..., is discrete, too. For each discrete value of 
the transverse momenta kj_ = (0, k y , k z ) T one then has 
four scattering states which we label <J> c ,kj_± and 3>h,kj_ ±j 



*h,kL,±(r) 



2e ±f M<0*-H*t,» sin{k z z) 
^/v x dW y 



^v x dW y 



(40) 



where k x {e) is the positive solution of 



; / \2 i 2 i 2 i 2 ^TfiS 

k x (E) = k F H — k y — k z -\ ^2~, 

and 



v x = hk x /m. 



(41) 



(42) 



(43) 



The scattering states labeled "+" represent quasiparticle 
states moving to the left; the states labeled "— " represent 
quasiparticle states moving to the right. All scattering 
states are normalized to unit flux. 

In the region x > 0, the scattering states differ from 
those given above because of the finite tunnel coupling to 
the superconductor. In particular, the scattering states 
acquire a finite weight inside the superconductor. In the 
tunneling limit r< 1, this weight is small and the ma- 
jority component of the scattering states inside the half 
metal remains well approximated by Eqs. ([4T)| and (|4"Tj) 
above. The exact expressions for the full scattering state 
in the region x > are cumbersome, and we refer to Ref. 
ll7L where the detailed expressions can be found. 

The "turning on" of the tunnel coupling to the super- 
conductor at x = gives rise to a small amount of normal 
reflection, but it does not cause Andreev reflection. We 
neglect this normal reflection at x = in the remainder 
of this section. 



B. Andreev reflection from a superconducting 
segment of length 8L 

The presence of spin-orbit coupling in the supercon- 
ductor gives rise to Andreev reflection at the HS inter- 
face, as we have seen in Sec. IIII1 In the second step 
of our calculation, we compute the effective Andreev re- 
flection amplitude for a superconducting segment of size 
< x < SL. We choose the length 6L of the supercon- 
ducting segment such that d,SL <C £. The inequality 
d <C £, together with translation symmetry in the y di- 
rection, ensure that the Andreev reflection amplitude is 
diagonal in k y and k z . The inequality 5L <C £ gives 
\k x (e) — k x {— e)\5L <C 1 for excitation energies up to 
Ao- This, in turn, leads to Andreev reflection amplitudes 
proportional to SL, which we write as p^ e (k±,£)5L and 
A=h(kj_, e)SL, for electron-to-hole and hole-to-electron re- 
flection, respectively. 

Calculating the Andreev amplitudes for the segment 
< x < 5L in first-order perturbation theory in the 
spin-orbit interaction gives 

/9 h e(kj_,e)(5i = -i($h,k L ,-|#WsoI$e,k u +) ) 
Peh(k±,e)SL = -£($e,kx,-|S%so|$h,k-L,+), ( 44 ) 
where SUso is the 4x4 matrix Hamiltonian representing 
the projection of the spin-orbit interaction Hamiltonian 
onto the segment < x < SL, 

l{P SL (x),H so }, (45) 



sn 



so 



%so — 



Hso 





-H, 



so 



(46) 



with P$l{x) = 1 for < x < SL and otherwise. Eval- 
uating the matrix element in the limit d <C SL <C £ then 
gives 



Phe(k_i_,e) = r hc (k||,e) 



(47) 
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FIG. 3: A lateral contact between a half-metallic film or wire 
and a superconductor. The contact has length L, with a seg- 
ment of length SL singled out. The main text describes how 
the effective Andreev reflection of the entire lateral contact 
is calculated, starting from the effective Andreev reflection 
amplitudes phcSL and p e h5L of the segment of length SL. As 
discussed in Sec. IVI1 inserting a normally reflecting boundary 
at the left end of the H film or wire gives rise to a Majorana 
bound state. 



where kii = (k x (0), k y , 0) T . Equation (|4"T)) has the sim- 
ple interpretation as the Andreev reflection amplitude 
for a single reflection at the HS interface, multiplied 
by the number of bounces at the HS interface per unit 
length . 17 ' 18 Similarly, one finds that 

p ch (k ± ,e)=r ch (kj | , £ )^-, (48) 

where kj, = (-fc x (0), k y , 0) T . 

In the same way, one also calculates Andreev reflection 
amplitudes Phc(kj_, e)'5L and Phe(kj_, e)'5L for quasipar- 
ticles incident on the segment < x < SL from the right. 
These are 



Phe(k_L,e)' = r he (kj|,e) 
p Gh (k_L,e)' = r c h(k||,e) 



2dk 
2dk 



(49) 
(50) 



Effective Andreev reflection amplitudes for the 
lateral contact 



In the final part of the calculation, we consider a super- 
conducting contact of finite length < x < L. In order 
to keep the notation simple, we omit the arguments k^ 
and e in the intermediate results. 

Upon comparing contacts of length L and L + SL, one 
finds, that 
rf e (L + 5L) = r f c ( L y(k,(e)-k^-e))SL 

-p' ch rf c (L) 2 ]SL + 0(5L) 2 .(5l) 



+ [Ph. 

Making use of the relations 



Ph 



= -(Peh)*, Peh = -(Phe)*, 



(52) 



which follow from quasiparticlc conservation or from the 
explicit solutions obtained above, and expanding k x (e) « 
fc(0) + e/Tiv x plus terms of order e 2 that are neglected in 
the Andreev approximation, one arrives at a nonlinear 
differential equation for rj^(L), 



ar hc 
dh 



2is 

hv x 



r hc + Phe ~ Phc( r hc) 2 - 



(53) 



Solving this equation with the boundary condition 
rgf (0) = gives 



rf c (L) 



p he sinQhoL 



Qhe cos QhcL - i(e/hv x ) sin Q^ e L 
and, similarly, 

eflF/r\ pch sin Q cil L 

ch Q e h cosQ hi + i(e/hv x ) sin Q e hL' 

where we abbreviated 



Qhe =V (s/hv x ) 2 - |phe| S 
Qeh =V {s/frv x ) 2 - |p ch | 5 
eS I 



(54) 



(55) 



(56) 



Upon taking the expression for r^(L) to first order 
in L, one verifies that one reproduces the starting point 
rfz(8L) = phoSL of the previous subsection. Upon taking 
r ho (L) to first order in pi lc , one finds 



l(L)= Phe [ L ' dxe 2Ux ' hv * 

JQ 



(57) 



which one obtains by applying first-order perturbation 
theory to the entire superconducting contact of length L 
at onceji 7 - Equation ([57]) represents the effect of a single 
Andreev reflection at the HS interface, with the phase 
factor accounting for the relative phase shift between the 
electron and the Andreev reflected hole for an Andreev 
reflection taking place at position x. One checks for con- 
sistency that r h ^(SL)/5L turns into Eq. (|4"T)) in the limit 
SL -> 0. 

The relevant limit for the lateral contact of Fig. [TJ) is 
the limit L — > oo. In this limit, the energy dependence 
of rj^J through the energy dependence of phe can be ne- 
glected in comparison to the explicit energy dependence 
in Eq. (|54[) . so that one may approximate phe(k_L,e) by 
Phe(k^,0) in Eq. ([54]). Defining 



e (k±) = hv x \p hc (k ± ,0)\ 
7 >hc(k||,0)|, 



one then finds that 



2d 



p hc (k_L,0) 



(58) 



7 i arcsin(e/eo 



|Phe(k±,0)| 

^hc(k||,0) .i arcs in( E / eo (kj_)) 



Ko(k||,0)| 



(59) 



if £ < £ (kj_). Hence, as long as e < £o(kj_), |rg*| = 1 
such that a lateral contact serves as an "ideal" contact 
between a superconductor and a half metal, allowing per- 
fect spin-flip Andreev reflection back into the half metal. 
For e > £o(kj_), rho(k^, e) is an oscillating function of the 
contact size L, whereas the magnitude of the Andreev re- 
flection amplitude is a decreasing function of energy, 

sin 2 (Q hcJ L) 



C(kL lE )|' 



(e/eo(k±)) 2 -cos 2 (Q he i) 
1 - y/l - MkO/e) 2 if L 



(60) 
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where the last line is obtained by averaging L over a 
period < Qh c L < 2ir and is proportional to the en- 
velope of l^hcf ( e ) I 2 ■ Th° energy £o(kj_) separating the 
regions of complete and partial Andreev reflection can 
be interpreted as a mode-dependent proximity-induced 
"minigap" in the half metallic film. 

An analogous calculation including multiple ARs can 
be done for the transmission through the contact. For 
electrons and holes incoming from the left we find the 
transmission amplitudes 



C(k±,e) 



Ch(k±,e) 



ifc x (0)L 



Qhe cosQ ho L - i(e/hv x ) sinQ hc L 



she 



-ik x {0)L 



Qeh cos QchL - i(e/hv x ) sin Q c hL 



(61) 



(62) 



The amplitudes for particles incoming from the right are 
related to Eq. ((61J- E2]) by t' c f = e vk*(o)L t eB and t'^f = 
e -i2^(0)i t eff For a long contactj |p he (kj_,e)|L > 1, and 
e < £o(kj_) the transmission amplitudes can be approxi- 
mated as 



functional form of this decrease depends on the details of 
the spin-orbit coupling. For voltages much larger than 
the induced "minigap" in the half metal, but still much 
smaller than Ao, i.e., hv-p ^\rhc\ / 2d <C eV <C Ao, we may 
take \r^{k ±7 e)\ 2 from Eq. §U§ and find 



G{V) 



N(t 2 



128 



eVd 



(66) 



This decay of the subgap conductance with the applied 
voltage is a marked difference with the case of the serial 
geometry, for which G is an increasing function of V . 



B. Josephson current 

The expression for the Josephson current / in a lateral 
SHS junction can be obtained from Eq. (|32[) upon setting 
r cc = 7*hh = and upon replacing r^e an d r e h by r?* and 

'eh > 



respectively, 



C(kL,e) « 2e tk ^ 0)L e- q ^ L e 



i arcsin 



(e/eo(k±)) 



(63) 



with ghe = iQhc- Thus, i ee , thh become exponentially 
suppressed consistent with a fully developed proximity 
effect. 



V. APPLICATIONS: SUBGAP CONDUCTANCE 
AND JOSEPHSON CURRENT 

The effective Andreev reflection amplitude can be 
used for a calculation of the subgap conductance G and 
the Josephson current / in a lateral HS or SHS junction 
in the same way as the Andreev reflection amplitude The 
is used in the case of a serial junction. 



A. Subgap conductance 

In the limit of a long ballistic contact and at the Fermi 
level e = 0, the effective Andreev reflection amplitude of 
Eq. (|59|) has modulus 1 for all transverse channels (la- 
beled by the integer n and the wavenumber k y ). Hence 
the zero bias conductance G(0) is 



G(0) = G m = ^-N, 
h 



where 



N 



k^pWyd 

47T 



(64) 



(65) 



is the number of propagating modes at the Fermi level in 
the half metal. 

Upon increasing V, the Andreev reflection probabili- 
ties and, hence, the conductance G decrease. The precise 



4ek B TW v d„ d ^ 



4tt 2 



n=0 



dk. 



_ n J k±<k f 



x In 



1 - »hf (k_L, iu n )r e £ (k ± , iu n )'e 



(67) 



where Lj is the junction length, see Fig. [5p, r£f (kj_,e) 
is the effective electron-to-hole Andreev amplitude of 
the right superconducting contact for quasiparticles inci- 
dent from the left, and ^^(k^e)' the effective hole-to- 
electron Andreev amplitude of the left superconducting 
contact for quasiparticles incident from the right. We 
have set the phase of the superconducting order param- 
eter for the right superconductor equal to zero. Upon 
setting e = iuj ni the effective Andreev reflection ampli- 
tudes r£^(k_L, s) and (kj_, e)' have a well-defined limit 
for the contact size L — > oo, which is given by Eq. (|59[) 
for all uj„. 

We first consider the limit when the minigap £o 3> -^Th 
is much larger than the Thouless energy £?xh = ^z/Lj of 
the junction (long-junction limit). For high temperatures 
ksTLj/hUF 3> 1, one finds 



el 



2irG m k B Te B " F . H 



27rfc R TL i 



sin < 



(68) 



This expression for / has the same sinusoidal phase de- 
pendence and exponential junction length dependence as 
in the serial geometry, but in the lateral contact geom- 
etry / is proportional to the much larger conductance 
G m , Eq. <J64J), instead of Eq. ([30]). In the limit of zero 
temperature, the Josephson current is given by 



2e 7Vi>f,h 
~3~ U 



4>, — 7r < 4> < ir . 



(69) 
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and I(4> + 27r) = I{<j>). Equation (|69|) is the known form 
of a supercurrent if the superconductor and the normal 
junction material are strongly coupled and the junction 
is disorder-free^ The phase dependence is sectionally 
linear (saw-tooth-like) and the critical current decreases 
with the length Lj of the junction. 

In the opposite limit of a short junction (eq <C -Erh), 
we distinguish three temperature regimes. For very high 
temperatures, k B T i?Th = hvp^/Lj, one obtains a 
sinusoidal phase dependence of the Josephson current 



I = 



3eN(T 2 



iTr 2 hd 2 k B T V 2irk B T 



Et\i 



4 Th sin < 



(70) 



For intermediate temperatures, £o <C k B T <C E^h, one 
finds 

(71) 

For T — 0, the trace over transverse modes could not be 
performed in closed form. However, the dependence on 
phase difference <j) can be found. One obtains 

1= ^TrK|r hc |}sin^, -tt < cf> < n , (72) 

and + 2-k) = I(4>). The dependence / oc sin((/>/2) 
of the zero-temperature Josephson current is reminis- 
cent of the "fractional Josephson effect", characteristic 
of Josephson junctions that have Majorana bound states 
at the superconductor interface d 37 ! 53 (see also the next 
section). 



VI. MAJORANA STATES 

Majorana bound states have been proposed as an el- 
ementary building block of a topological quantum com- 
puter since they are an example of an excitation with 
non-Abelian statistics^ Majorana bound states exist as 
the fundamental excitations of a candidate state for the 
v = 5/2 quantum Hall effect ) 54 ' 55 in vortices in super- 
conductors with a spinless p-wave pairing symmetry^ 5 - 5 -"— 
or in vortices of s-wave superconductors in contact to 
a topological insulator— or a standard two-dimensional 
electron gas in a large magnetic field and with strong 
spin-orbit coupling^ 

Very recently, it was suggested that Majorana bound 
states can be found at the ends of semiconducting 
quantum wires with strong spin-orbit coupling and 
a strong magnetic field, in contact to an s-wave 
superconductor j 40 i 41 i 43 In these proposals, the role of the 
magnetic field is to create a gap for spin excitations, so 
that the wires become effectively half metallic. We now 
show that Majorana bound states also occur in the sys- 
tem considered here: a half-metallic quantum wire in con- 
tact to a superconductor with spin-orbit coupling. This 



enables us to make contact between the (experimentally 
observed) triplet proximity effect and the so far purely 
theoretical search for avenues to topological quantum 
computation. Our approach has the additional benefit 
of providing a fully microscopic description of the p-wave 
proximity state, in contrast to the existing studies of this 
effect in semiconducting wires with strong spin-orbit cou- 
pling, which rely on an effective description using an in- 
duced pairing potential in the semiconducting wirei 40 ' 41 
First, we show that a Majorana state exists at the end 
of a ballistic half-infinite half metallic quantum wire lat- 
erally coupled to a superconductor. We consider the ge- 
ometry shown in Fig. [3l The difference with the calcu- 
lation of Sec. [V] is that here the half metal is a wire, 
not a film. We therefore have to use hard-wall bound- 
ary conditions in the y-dircction, not periodic boundary 
conditions as in Sec. [Vj With hard- wall boundary condi- 
tions, the Andreev reflection amplitudes phc and p e h per 
unit length have to be replaced by amplitudes p\ lc and 
p c h, which are defined as 

Phc(k y ,k z ,e) = - [phc{k y ,k z ,e) + Phc{~k y ,k z ,s)} ,(73) 

Pch(k y ,k z ,e) = - \fieh(k y ,k z ,s) +p<ih(-ky, k z ,e)} .(74) 

Since phe and phc are odd in k y , the components of the 
SOI coupling to k y drop out. Apart from the replacement 
Phc Phc and Pch — > Pch; the results of Sec. IIVI continue 
to hold for the present case. 

At the left end of the half-metallic wire, quasiparti- 
cles undergo normal reflection with amplitude sh(£) and 
sh(~ e)* for electron-like and hole-like quasiparticles, re- 
spectively. With the Andreev reflection amplitudes of 
Eqs. (|59|) we then find a non-degenerate bound state at 
e = with (unnormalized) wavefunction 

*(r)=i[r£]"* [s H (0)$ c , ki , + (r) + $ o , ki ,_(r)] 

+ * [r^Y [$ h) kx,-(r) + s H (0)*$ h , ki , + (r)] , (75) 

where the scattering states $ ,kj_,± and $h,kj_,± are ob- 
tained from those given in Eqs. (|40| and (|4Tj) above, but 
with the replacement e lkyV — > sm(k y y) because of the 
hard-wall boundary conditions. The distance to the next 
bound states is of the order of the minigap £o or the 
level spacing in the normal segment extending from the 
superconductor, whichever is smaller. 

The bound state ([75)1 is identified as a Majorana bound 
state because it is invariant under electron-hole conjuga- 
tion, i.e., Ti<i>* = $ where t\ is the 1st Pauli matrix in 
electron- hole space. Alternatively, with ^1,^4 being 
electron and hole creation operators, $ corresponds to 
the field operator 



1=1 dx[u t (x)>$ t (x) + u t (x)^j(x)] 



(76) 



with u-f — $1 and = $3 given by the electron- and 
hole spin-up component of <I>, respectively. This operator 
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satisfies the condition 



7 = 7 f , 



(77) 



which is the defining characteristic of a Majorana state. 
Being a Majorana bound state, $ is stable against per- 
turbations because, by particle-hole symmetry, a pertur- 
bation that moves i> to some finite energy £ / must 
generate a pair of states at ±e. Since $ is a single state 
this is not possible. 

We note that there is one Majorana mode per trans- 
verse mode in the half-metallic wire. Disorder, which 
is not included here, will lead to interactions between 
these modes, which will cause Majorana modes to pair- 
wise combine into standard fermionic excitations. If the 
number of transverse modes is odd, a single Majorana 
mode is guaranteed to remain present at the end of the 
half-metallic wire42 

If the half-metallic quantum wire has a finite length 
L, the Majorana bound states at the two ends will in- 
teract, so that the excitation acquires a finite energy, 
exponentially small in the length of the wire. This finite 
excitation energy can be calculated from the full scatter- 
ing matrix S(e) of the lateral H-S contact, calculated in 
Sec. [iVl and the reflection amplitudes sn(e) and s' H (e) 
at the left and rights ends of the half-metallic wire. The 
energy spectrum is found from the condition 



where 



ShOO = 



S(s) 



det(l-5 H (e)5(e)) = 0, 



s H (e) 

s H (-e)* 

s' a (e) 

s' K (-e) 

( (e) tf c (e) \ 



(78) 



Ch(e) 



£e(e) rf( £ ) 
V tll(e) 



In the limit of a long contact, |/5h |£ 3> 1, we then find 



= ±2\p he (e)\nv x (e)e-\"^ L \ sin(fc x (e)L)| 



(79) 



where we have set sh = Sjj = — Thus, the energy split- 
ting decreases exponentially with the contact length (be- 
sides accidental degeneracies for integer k x (e = 0)L/2n). 

It is instructive to compare our calculation with 
the model of a spinless one-dimensional p-wave 
superconductor,— which has been used as a phenomeno- 
logical model description of the induced superconductiv- 
ity in a semiconductor wire with a strong magnetic field 
and spin-orbit coupling . 39 ! 42 This model has the Hamil- 
tonian 



H 



P 

2m 



ro + A'pn, 



(80) 



where To is the 2x2 unit matrix in electron- hole space and 
A' the effective p-wavc superconducting order parameter. 
Comparing with our calculation, and specializing to a 
quantum wire with one quantized mode only, for which 
k z = Tr/d, we identify 



lA'l 



irhT(k±) 



4d 2 (k 2 



F,S 



1.2 _ U2\ 



"1 



o 2 



(81) 



where T(kj_) is the transparency of the interface at the 
relevant (lowest) transverse mode. 



VII. DISCUSSION AND CONCLUSION 

In this article, we have shown that spin-orbit interac- 
tion in a singlet superconductor gives rise to a triplet 
proximity effect if the superconductor (S) is coupled to 
a half-metallic ferromagnet (H). We have calculated the 
conductance of a HS junction and the Josephson current 
of a SHS junction in both a serial geometry and in a lat- 
eral contact geometry. Because of the coherent effect of 
multiple Andrcev reflections, the effective Andrcev am- 
plitudes for a lateral contact geometry are significantly 
enhanced in comparison to those at a serial geometry. In 
particular, multiple Andreev reflections at the interface 
between a clean (disorder-free) half-metallic film or wire 
and a superconductor can lead to a fully developed triplet 
proximity effect in the half metal, with an Andreev re- 
flection amplitude of unit magnitude. 

The results found here have been derived under the as- 
sumption of a ballistic system, i.e., without taking into 
account disorder scattering in the half-metal or in the 
superconductor. For the single-reflection amplitude The 
in Sec. Ill Bl this does not strongly restrict the validity 
of the result since the Andreev reflection amplitude is 
a microscopic property of the interface: is deter- 
mined by matching the eigenstates on length scales of 
the Fermi wave length in H and S and the wave func- 
tion decay length in S. If the disorder is weak, such that 
the mean free path I exceeds these microscopic length 
scales, The will be unchanged by the presence of disorder. 
In this way, the microscopic Andrcev reflection ampli- 
tudes calculated in Sec. Ill Bl mav also serve as a starting 
point for studies of the conductance of a disordered HS 
junction or the Josephson current in a disordered S-H-S 
junctions. (For example, the Josephson current through 
a disordered or a chaotic Josephson junction can be found 
by combining the reflection amplitudes of the clean su- 
perconductor interface with the normal-state scattering 
matrix of the junction, see, e.g. Refs. I50ll5lh 

On the other hand, quantities that rely on free (phase- 
coherent) propagation in Hmay change qualitatively in 
the presence of disorder. Specifically, the effective An- 
dreev reflection amplitude r^f of the lateral contact has 
been obtained by phase-coherently summing single re- 
flection amplitudes. At the Fermi energy, these multi- 
ple Andrcev reflections add constructively because the 
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momentum k|| parallel to the interface is conserved and 
amplitudes of subsequent reflections have the same sign. 
However, scattering from impurities under the contact 
will lead to a summation over single amplitudes with 
different incident angles. Since the Andreev reflection 
amplitudes rho and r h are odd in k|| , this sum may no 
longer be constructive. Thus, the result for r^f is valid 
for an ideal, disorder-free lateral contact only. Since the 
reflection properties of a lateral junction saturate if the 
junction length L > d/|rh e |, where d is the thickness of 
the half-metallic film or wire, disorder is not expected to 
significantly alter our results as long as the elastic mean 
free path I ^> c?/|r ho |. 

As a particularly timely application of our calculation, 
we connect the Andreev reflection amplitudes and 
r c h calculated here to the existence of Majorana bound 
states at the ends of a ballistic half-metallic quantum wire 
in (lateral) contact to a superconductor with spin-orbit 
coupling. This proposal for the construction of Majorana 
bound states is a variation of a recent proposal that such 
Majorana bound states exist at the ends of a semicon- 
ducting wire in contact to a superconductor, where the 
semiconductor has strong spin-orbit coupling and the sys- 
tem is placed in a large Zeeman field, such that a the wire 
becomes effectively half metallic ! 40 ' 41 In our construction, 
the Zeeman field is replaced by the exchange field in the 
half metal, and the spin-orbit coupling is not located in 
the wire, but in the superconductor. It thus avoids the 
necessity of a (fine-tuned) applied magnetic field, which 
could negatively interfere with the superconducting or- 
der. 

We thank B. Bcri, J. Kupfcrschmidt, D. Manske, O. 
Starykh, A. Schnydcr, F. von Oppcn, and F. Wilken for 
discussions. This work was supported by the Alexander 
von Humboldt foundation. 



Appendix A: Perturbation theory in Hso 

An alternative calculation of the Andreev reflection 
amplitudes r^ c and r G h makes use of perturbation theory 
in the spin-orbit interaction iTso- For this calculation, 
finite lateral dimensions W x x W y are assumed, with pe- 
riodic boundary conditions in the x and y directions. 

As before, we consider wavefunctions proportional to 
e ik m x+ik y y _ j n ^ ne absence f spin-orbit coupling, there 
are two linearly independent solutions of the Bogoliubov- 
de Gennes equation for each pair k x ,k y . The first of these 
is "electron-like" , and of the general form 



*ek„ (r) 



Cef e 



ik z (e)z I / 



+ c ct e 



-ik z {e)z 



y/v z (e)W x W y 
x e^ x+ifc " 2/ (l,0,0,0) T 
+ c hi e tk ^ x+lk y v+K ^-^ z (0, 0, 0, 1) T (Al) 



for z < and 

#ek„(r) = d\e lk * x+lk * v+tq+z (1,0,0, e -^-"?) T 

+ ^e^+ikyy+iq-^ ^ 0) 0) e -«/>+»)) T (A2) 

for z > 0, where 



s\/k^ s -k x -k y + 2ism x /Al-£ 2 . (A3) 



The boundary conditions (|22[) with Qj >z — give four 
equations for the five coefficients c e ,f, c' e p Ch.j,, (if-, and 
dp so that one coefficient can be chosen freely. Choos- 
ing Cerf = 1 one obtains the "retarded scattering state" 
Vlf^k , while choosing cf e ^ = 1 one obtains the "advanced 

scattering state" ^f^^. . 

e,K|| 

The second scattering state is "hole-like" and has the 
general form 



*hk„(r) 



y/v z (-e)W x W v 
x e lk * x+lk y y (0,0, 1,0) T 
+ CBj> e tt - x - H *»* + "« (e) * (0, 1, 0, 0) T (A4) 



for z < and 



#hk u (r) = d! i e ik ' x+ik y y+iq+z (0,l,-e-^ t?) ,0) T (A5) 



for z > 0. The boundary conditions (|22|) with = 
give four equations for the five coefficients c c ,±, Ch,f, c' h p 
d±, and d^> . (See Ref. [13 for details.) Choosing c ,t = 1 
one obtains the "retarded scattering state" ^^ky > wnnc 
choosing c' e + = 1 one obtains the "advanced scattering 
state" . 

e,k|| 

In the Born approximation, the Andreev reflection am- 
plitudes to first order in the spin-orbit interaction are 
then found as the matrix element 

r ho (k||, £ ) = -*« k|! |Hso|*e R k ll >, (A6) 
r ch (k||, £ ) = -*(* c A k|| |H S o|^ k|| ), (A7) 

between retarded and advanced scattering states, where 
%SO is given by Eq. (|46| . 

Inserting the explicit expressions for the scattering 
states into Eqs. (|A6|) and (|A6|) then gives the Andreev 
reflection amplitudes of Eqs. (f2"3")l and (|M|) . 

Note that only the wavefunctions in the superconduct- 
ing region z > enter into the calculation of the Andreev 
reflection amplitudes. The observation that the SOI pro- 
portional to f2i iZ or Q2,z does not give rise to Andreev 
reflection to first order in the SOI then follows from the 
observation that the matrix elements (|A6j) and (|A7jl van- 
ish for arbitrary coefficients d^, dp d^, and d'^ if -ffso 
contains f2i iZ or SXi.z only. 
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Appendix B: Triplet pairings in S 



and 



In this appendix we give the results for the Andreev 
reflection amplitudes in the presence of finite triplet pair- 
ings in S which are linear in momentum. The triplet 
proximity effect in a half-metal in contact to a spin-triplet 
superconductor was also considered by Linder et al*^ We 
consider a pairing of the form 



A i (p) = J2d i 



Pi 



(Bl) 



where dij are the components of the triplet order pa- 
rameter which is taken as a small correction to the spin 
singlet s-wavc pairing Aoe 1 ^ (see Sec. [TTJ) . In a similar 
calculation as in the main text we find the Andreev re- 
flection amplitude induced by Aj(p) to first order in 



Aj(p). The amplitude r^ e 



' he 



r h ' e contains contri- 



butions r 



he ' ho 



that are even and odd in momentum, 



respectively We find 



he 



gr^n ^ + ^) V*f,s - k l.H «n 2 9 (B2) 



2(A2- £ 2) 



' he 



3 

3=1 



ir(0) 



d *n + id h 



A§ (d* l3 + id* 2j ) 



(d lj +id 2j ) 



2 (Ag-e*)* 
for the electron-to-hole conversion amplitudes and 
^eh(k||,£) = K (-k||,-e)]* 



(B3) 



(B4) 



for the opposite process. Triplet pairings can be included 
in the calculation of the conductance and the Josephson 
current by the substitution 



'he 



' he ■ 



(B5) 



The enhancement due to multiple Andreev reflections 
found in the second part of the article for the lateral 
geometry does not rely on the details of The and is not 
changed by the presence of the Aj. For a detailed dis- 
cussion of amplitudes that are odd or even in frequency 
we refer to Ref. 
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